Shubnikov - de Haas effect in the quantum vortex liquid state of the organic 
superconductor k-(BEDT-TTF) 2 Cu(NCS) 2 
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We report the Shubnikov-de Haas (SdH) oscillations observed in the vortex liquid state of the quasi 
two dimensional organic superconductor k-(BEDT-TTF)2Cu(NCS)2- The SdH oscillations can be 
observed down to about 5 T at 0.5 K, where the flux flow resistivity becomes as small as about 30 % 
of the normal state value. Below the upper critical field H C 2 of about 7 T, the additional damping of 
the SdH oscillation amplitude appears, as well as that of the de Haas- van Alphen (dHvA) oscillations, 
with respect to the normal state one which is described with the standard Lifshitz-Kosevich formula. 
The magnitude of the additional damping near H C 2 is the same with that observed in the dHvA 
oscillations and well explained by the theoretical predictions in consideration of fluctuations in the 
thermal vortex liquid state. In the quantum fluctuation region at lower temperature, however, only 
SdH effect shows the stronger damping than that of the dHvA oscillations. The different magnetic 
field dependence of the additional damping of the oscillation amplitude between the SdH and dHvA 
effects is discussed in connection with the effect of the transport current on the short-range order 
of vortices in the quantum vortex slush state reported at the same temperature and magnetic field 
region. 

PACS numbers: 74.70.Kn, 71.18.+y, 74.40.+k 



I. INTRODUCTION 

After the report of the magnetic quantum oscillations 
in the superconducting state of 2if-NbSe 2 more than 
quarter century agojj the oscillations of the magnetiza- 
tion, de Haas - van Alphen (dHvA) effect, in the vortex 
state seems to be confirmed experimentally for a [Variety 
of the type-II superconductors in the last decade.0 Com- 
mon experimental results of the dHvA oscillations in the 
superconducting state are that the additional damping of 
the oscillation amplitude appears below the upper criti- 
cal field H C 2 with respect to the normal state damping. 
The additional damping has been discussed in several 
ways theoretically.El The questions in those discussions 
are summarized as follows ;u (1) if the superconducting 
gap Ao exists in the vortex state just below i/ c2 , it would 
have drastically damped the oscillations at low temper- 
atures by a factor of exp (— A /£;bT), (2) the inhomo- 
geneous field distribution due to the flux lattice would 
broaden the Landau levels, and (3) the inhomogeneity 
in the superconducting order parameter associated with 
the vortex lattice leads to inhomogeneous broadening of 
the Landau levels in the quasiparticle spectrum near the 
Fermi surface. Thus the oscillations include basically rich 
information on the quasiparticle in magnetic fields and 
also the vortex matter properties .□ 

Shubnikov-de Haas (SdH) oscillations in the supercon- 
ducting state is very difficult to be observed because the 
finite resistivity is needed inevitably. Then it can ap- 
pear in the quite limited field-temperature region where 
the long-range translational order of the vortices is lost. 
Quasi two dimensional (Q2D) organic superconductors 
are good candidate for the observation of SdH effect in 
the vortex state because Affge fluctuations induce the 
wide vortex liquid region.BQ At low temperature below 



1 K, the quantum vortex liquid (QVL) is realized in k- 
(BEDT-TTF) 2 Cu(NCS) 2 (T c ~ 10 K) due to the large 
quantum fluctuation instead of the thermal one.l3tl In the 
QVL region, the finite resistivity is expected to remain 
even below H C 2(T ~ 0). Recently the transport proper- 
ties in the QVL region have been examined in detail, and 
the finite resistivity has been confirmed!] 

In this paper, we report the SdH oscillations observed 
in the QVL region of the Q2D organic superconductor 
k-(BEDT-TTF) 2 Cu(NCS) 2 . An additional damping of 
the oscillation amplitude appears on both the SdH and 
dHvA effects around H C 2, which may come from the su- 
perconducting fluctuation. At lower temperature in the 
QVL region, however, the stronger damping is observed 
only on the SdH effect. The different magnetic field de- 
pendence of the additional damping of the oscillation am- 
plitude between the SdH and dHvA effects is discussed 
in connection with the effect of the transport current on 
the short-range order of vortices in the quantum vortex 
slush statea reported at the same temperature and mag- 
netic field region. 



II. EXPERIMENT 

High quality single crystals of k-(BEDT- 
TTF) 2 Cu(NCS) 2 were grown by an electrochemical 
oxidation method. The magnetic torque measurements 
were performed by using precision capacitance torqueme- 
tcr. The in-plane and the out of plane resistivities were 
measured along the b and a* axes, respectively, by means 
of a conventional ac or dc four terminal method. The 
electrical terminals were made of evaporated gold films, 
and gold wires (10 fiva) were glued onto the films with 
gold or silver paint. The contact resistance was about 10 
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FIG. 1: Magnetic torque curves of Fi lHe T §S» 1 i i p4e a i#l in «- 
(BEDT-TTF) 2 Cu(NCS) 2 at 0.52 K and in the field direction 
tilted by 9 =6.1° from the a* axis, (a) the dHvA oscillations 
around 7 T in the expanded scale, (b) the magnitude of the 
torque hysteresis. 



il for each contact at room temperature, but it became 
less than 1 Q at low temperature where the experiments 
were carried out. The torquemeter and the samples for 
the resistivity measurements were fixed to the single axis 
rotation holder which can change the sample direction 
with respect to the magnetic field with the accuracy of 
0.05 degree. The holder with the samples was cooled 
slowly from room temperature to 4.2 K in 48 hours and 
specially slow cooling rate was used between 50 and 75 K 
in order to avoid the disorder of the— terminal ethylene 
group of the BEDT-TTF molecules.EEl The holder was 
directly immersed in liquid 3 He of the refrigerator which 
was combined with a 15 T superconducting magnet at 
the High Magnetic Field Laboratory for Superconduct- 
ing Materials (HFLSM), IMR, Tohoku University. The 
results presented in this paper were obtained on three 
samples #1 for the magnetic torque, and #2 and #3 for 
the resistivity measurements from different batches. We 
found that other one sample for the magnetic torque and 
one sample for the resistivity gave qualitatively similar 
results which were not presented in this paper. 



III. RESULTS AND DISCUSSION 

Figure 1 shows the magnetic torque curves of the sam- 
ple #1 at T — 0.52 K. The overall features are the same 
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FIG. 2: Magnetic field dependence of the resistivity pb and 
p a along the b (in the Q2D plane) and the a* F ixes il foui; l of the 
plane), respectively, in the sample #2. The upper left inset 
shows the first Brillouin zone and the Fermi surface. The 
lower right inset indicates the SdH oscillations on p a plotted 
as Aa/a n - 

with the previous report ;@ the irreversible and reversible 
regions are separated at H m ~ 3 T (Fig. 1(b)), and 
the dHvA oscillations with one fundamental frequency of 
F a = 599 ± 2 T are observed in both the normal and 
superconducting states. The quality of the sample used 
in this study seems to be better than that in the previous 
world] judging from the large amplitude of the dHvA os- 
cillations. The reversible magnetic torque region (Hi„ ~ 
3 T < H < H c2 ~ 7 T) below 1 K is expected to be the 
QVL region. The finite resistivity appears in the QVL 
region even at T ~ 0. The detail of the transport prop- 
erties in the QVL region has been already reported^ It 
is noted only here that a weak non-linear behavior of the 
resistivity is found in the QVL region. Such non-linearity 
is not observed in the thermal vortex liquid (TVL) region 
above 1 K. The concept of the quantum vortex slush has 
been proposed for the non-linear behavior below 1 K. The 
vortex slush with only the short-range <yi/?r. of vortices 
has been found in the high-T c oxides Mll3£j The effect 
of the quantum vortex slush and the TVL states on the 
dHvA and SdH effects will be discussed latter. 

Figure 2 shows the magnetic field dependence of the 
resistivity in the sample #2 at T = 0.52 K. The in-plane 
(pb) and the out of plane (p a ) resistivities are measured 
in one single crystal along the b and a*-axes, respec- 
tively. The magnitude of the magnetic field dependence 
of the resistivity and the resistivity onset corresponding 
to i/i rr are not influenced by the applied current density 
in both configurations. The flux flow resistivity in the 
QVL region, however, changes a little with the current 
density because of the weak non-linear resistance!!! The 
two curves in this figure are measured using the current 
density (current) of J = 0.16 A/cm 2 (/ = 100 //A) for p?, 
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and 1.6 mA/cm 2 (10 pA) for p a , respectively. The SdH 
oscillations with the frequency of 599 ± 2 T are clearly 
observed. The oscillations come from the a-orbit cen- 
tered at the Z-point of the first Brillouin zone depicted 
in the upper left inset. In higher magnetic fields, the 
magnetic breakdown orbit (3, consisting of the a and 
opened 7 orbits, has- been observed in both SdH and 
dHvA oscillationsoiij We, however, restrict ourselves to 
the single band model for the following analysis and dis- 
cussion on the SdH and dHvA effects because the mag- 
netic field used in the presept study is smaller than the 
magnetic breakdown fieldJlirta The amplitude of the SdH 
oscillation on p a is much larger than that on pf,, although 
the magnitude of the magnetic field dependence of the re- 
sistivity is almost the same with two configurations. The 
reason is not known at present but the similar tendency 
is commonly seen in the Q2D organic conductors.Ej 

In order to see the SdH oscillations in low magnetic 
field region, Aa/a n is shown in the lower right inset of 
Fig. 2. Here, Act is the oscillatory part of the conduc- 
tivity obtained by subtracting the non-oscillatory part of 
the conductivity, and cr n (~ p~ 1 ) is the normal conduc- 
tivity which is the same with the non-oscillatory part of 
the conductivity in the normal state. In the vortex state 
below about 7 T, SdH oscillations come from the nor- 
mal (quasi particle) component of the total conductivity 
which includes additional non-equilibrium superconduct- 
ing component due to the vortex pinning. The normal 
(quasi particle) part of the conductivity roughly corre- 
sponds to the flux flow conductivity. Then the normal 
resistivity p n {— "n" 1 ) in the vortex state is assumed to 
be extrapolated smoothly from the normal state into the 
vortex state toward p — at H = 0. The dotted line 
in Fig. 2 shows the normal (quasi particle) resistivity 
p n (~ cr^ 1 ) in the normal (vortex) state. The line is ob- 
tained as almost linear in the vortex state and passing 
through the middle of the SdH oscillations in the normal 
state. In intermediate region around 7 T, two lines are 
connected smoothly by the simple second order polyno- 
mial function. The SdH oscillation is persisting down to 
about 5 T where the resistivity is about 30 % of the nor- 
mal state value (~ 1.5 ficm at 7 T). In the other sample 
#3 the oscillations can be seen in lower magnetic field of 
about 4.5 T, which are shown latter in Fig. 4. This is first 
unambiguous observation of the SdH oscillations on the 
flux flow resistance in the well characterized supercon- 
ductor. Similar SdH effect in the superconducting state 
has been reported in the /3"-type BEDT-TTF organic 
superconductor .Lj But it did not show the additional am- 
plitude damping in both the SdH and dHvA oscillations 
in the superconducting state. It may be necessary to cop=. 
sider the smaller H C 2 value reported by another groupLa 
than the value expected in Ref. |l7|. 

The additional damping of the SdH and dHvA oscil- 
lation amplitude in the vortex state is demonstrated in 
Fig. 3. The oscillation amplitude At and Aa/a n for the 
dHvA and SdH effects is shown by the circles in Figs. 3(a) 
and 3(b), respectively. In the normal state above about 
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Fig.3 Sasaki et al. 

3: Magnetic field dependence of the oscillation ampli- 



tude of (a) the dHvA and (b) SdH effects, (c) The additional 
damping factor in the-pjrtex state /fee; Mr-by Maniv et alB 
and IC by Ito et al.ci and Clayton et a/.LJ The solid and 
broken curves in (a) and (b) are the field dependence of the 
oscillation amplitude Af^ and Alk based on the LK formula 
for the 3D and 2D cases, respectively. 



7 T, both oscillation amplitude are well described by the 
standard Lifshitz-Kosevich (LK) formula.EJ The ampli- 
tude A l k of the first harmonics of the oscillations in the 
single band is given by Alk <x TH u RtRdRs, where the 
temperature factor Rt = (Am c T / H) / sinh(Am c T '/ H) , 
the Dingle factor Rd = exp(— Xm^T^/H), and the spin 
factor Rs = cos(irgmi,/2mo). Here, A = 2ir 2 ckB /eh = 
14.69 T/K, m&, m c and tuq are the band, cyclotron 
effective and free electron masses, and Tr> is the Din- 
gle temperature related to the scattering rate tq with 
To = h/2nkBTQ. The power n of H depends on the mea- 
sured quantity and the dimensionality. For the SdH and 
torque-dHvA amplitude (Aer/<7 n and At) in the 3D (2D) 
case, n is 3/2 (1). The fitting in both cases are very good 
with m c = 3.5m and T D = 0.26± 0.05 K (0.17±0.03 K) 
for dHvA and 0.24 ± 0.05 K (0.16 ± 0.03 K) for SdH ef- 
fects in the 3D (2D) caseril The 2D formula, however, 
is adopted for the latter analysis. Because the quantized 
Landau level spacing in the present magnetic field region 
is considered to be fairly larger than the interlayes, trans- 
fer integral of this Q2D organic superconductor £3 

Below about 7T the amplitude starts to deviate 
smoothly downward from Alk in both dHvA and SdH 
effects. This indicates the additional amplitude damping 
in the vortex state, which has been reported so far on 
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the dHvA effect in K-(BEDT-TTF) 2 Cu(NCS) 2 .Qfia£ 
The origin of such smooth damping has been dis- 
cussed on the basis of the model of the quasiparti- 
cle scattering by the random vortex lattice with the 
large superconducting (vortex) fluctuation around t 
mean field fl^ffi In the approach by Maniv et al. 
the additional damping term Rsc as being multi- 
plied to A LK is R sc = exp(-7r 3 / 2 (|A| 2 )/n]/ 2 ), where 
np = Ep/htu c , tu c = eH/m c c, Ep the Fermi en- 
ergy, and (|A| 2 ) is the mean square of the supercon- 
ducting order parameter averaged over space coordi- 
nates. The tilde above A indicates that energy is mea- 
sured in units of Tilo c . This expression is very 
to the Maki Stephen type mean field approach] 2 -]! 2 ^ 
In ra simple analytic expression derived by Maniv et 
al.fi <|A| 2 ) ~ (q//3)[l + e xp(-x 2 )/2x exp(-y 2 )rfy], 

where a = (l/ 2fwj c ) In y/H c2 /H, /5L== 1.38 / n F (huj c ) 3 , 
and. x — a/y / 2f3kBT. Ito et a/.E£l and Clayton et 
al.ci have used an approximated interpolation formula 
for the mean square of the gap function, (|A| 2 ) = 
^[A(0) 2 (1 - H/H c2 )/2f + a(T) 2 +A(0) 2 (l-i// J ff c2 )/2, 
where a(T) is a temperature dependent parameter to 
scale the fluctuations. Both calculations of Rsc are 
shown in Fig. 3(c) by using the same mean field H C 2 = 
4.8 T. The results on i?sc are in fairly good agreement 
with each other. 

The dotted curves in Figs. 3(a) and 3(b) show the ex- 
pected magnetic field dependence of the oscillation am- 
plitude taking Rsc into account. The dHvA effect is 
well represented in .these fluctuation approaches as has 
been reportedEcH'Ei In the SdH effect, the experimental 
results follow well A^Rsc near iJ c2 ~ 7 T as well as 
the dHvA effect. A stronger damping, however, appears 
below about 6.5 T in only the SdH effect. 

We discuss the stronger damping in the vortex state 
observed only in the SdH effect in connection with the 
quantum vortex slush state at low temperature. In the 
QVL region the finite resistance appears in the vortex liq- 
uid state between the melting H m or irreversible field H UI 
and H C 2 even at T ~ 0. In the case of less or no quan- 
tum fluctuation, the zero resistance should appear just 
below H C 2 and H m (T) or H- m (T) coincides with H C 2(T) 
at T = 0. The former QVL region has been actually 
found in k-(BEDT-TTF) 2 Cu(NCS) 2 as a demonstration 
of the importance of the quantum fluctuations in this 
material .Era Recently two wortex liquid regions have been 
found at low temperature n The low resistivity state with 
non-linear current response below about 1 K has been dis- 
tinguished from the high resistivity state at higher tem- 
perature. A steep drop of the resistivity around Tl ~ 
1 K separates the vortex liquid state into these two re- 
gions. The short-range order of vortices has been ex- 
pected to exist in the former low resistivity state referred 
to as the quantum vortex slush state. Because these fea- 
tures in the low resistivity state are phenomenologically 
similar to the observations explained by the the vortex 
slush concept with the short-range order of vortices in 
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FIG. 4: SdH oscillations in the sample #3. Fi fnlef k sflows the 
overall magnetic field dependence of the resistivity p a at T = 
0.43 K. The dotted line indicates the normal (quasi particle) 
resistivity p n ~ a" 1 . 



the high-T c oxidesBBEl The latter high resistivity state 
has been considered as the thermal vortex liquid (TVL) 
state where no translational long-range order of vortices 
is formed. 

Since the additional damping of the oscillation ampli- 
tude in the vortex state is explained to be sensitive to 
the local phase modulation of the superconducting order 
parameter|] some kind of perturbation on the vortices, 
for example, applying the transport current in the SdH 
measurements, arc expected to influence the phase coher- 
ence. It is noted that both dHvA and SdH oscillations 
in this experiment are measured in the quantum vortex 
slush region. In the case of such measurements in the 
TVL state, the additional damping in the vortex state 
is expected to be the same in dHvA and SdH effects be- 
cause TVL does not have any order of vortices. Thus 
neither the presence of the transport current in the mea- 
surements of the SdH effect nor the absence of current 
in the dHvA effect alters any phase coherence. On the 
other hand, SdH effect in the quantum vortex slush state 
may be affected by moving vortices due to applying the 
transport current. It may disturb such coherency of the 
vortices and quasiparticles. Non-observation of stronger 
oscillation amplitude damping in the dHvA effect with 
basically no transport current demonstrates the influence 
of the moving vortices by the current in the quantum vor- 
tex slush state. 

Besides applying transport current, the effect of mov- 
ing vortices on the damping of the dHvA oscillations has 
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FIG. 5: SdH oscillation amplitude scaled By'Aj^ to'the sam- 
plers at 0.43, 0.76, and 0.88 K. Upper left inset shows thedbw 
temperature part of the schematic vortex phase diagram.H N 
is the normal, S the vortex solid, VS and TVL the quantum 
vortex slush and thermal vortex liquid states, respectively. 
Lower inset demopatrates the additional damping factor pro- 
posed by Ito et alxQ and Clayton et alX3 Each curve from top 
to bottom is calculated with a = 0.10, 0.15, 0.20, 0.25, and 
0.30 (meV) 2 . 



been seen in the superconducting state of 2.ff-NbSe2.0 
The damping of the oscillations depended on the history 
of reciprocal partial sweeps of the external magnetic field 
in the hysteretic region where the pinning strength was 
changing. For explaining the observations, Maniv et at 
have suggested that the motion of vortices, depending 
on the pinning strength, influences the magnitude of the 
damping of the oscillations 

In order to see the difference of the damping in the 
quantum vortex slush and TVL regions, we tried to com- 
pare the SdH oscillations at higher temperatures. Figure 
4 shows the SdH oscillations at higher temperatures in 
the sample #3. The SdH oscillations in the vortex state 
are clearly seen down to about 4.5 T at 0.43 K. The sam- 
ple #3 seems to have better quality to see the SdH os- 
cillations in lower magnetic field than that in the sample 
#2 presented in Fig. 1. 

The magnetic field dependence of the observed SdH os- 
cillation amplitude scaled by A^ in the normal state is 
shown in Fig. 5. The value of (Aa/a n )/Aj^ corresponds 
to the additional damping in the vortex state. The addi- 
tional damping starts smoothly around H C 2 ~ 7 T. The 
magnitude of the additional damping becomes larger in 
higher temperatures near H C 2- It is naturally understood 
by the fluctuations. The lower inset shows calculations 
of the additional damping factor i?sc proposed by Ito et 



alM and Clayton et aM The fluctuation parameter a is 
expected to take larger value in higher temperature. This 
has been actually confirmed in the dHjiA experiments 
on the same material by Clayton et alt3 The a value 
has changed almost continuously from 0.13 at 0.03 K to 
0.32 (meV) 2 at 0.44 K. It is noted that the R S c curves 
calculated with a or corresponding temperature are con- 
tinuously shifting to smaller value, and do not cross each 
other. The present SdH oscillation damping near H C 2 is 
in agreement with the theoretical calculations. But the 
oscillation amplitude at 0.43 K shows stronger additional 
damping in low magnetic field than others measured at 
high temperature. This stronger damping observed at 
0.43 K is considered to be due to the moving vortices in 
the quantum vortex slush state as discussed above. 

It is expected that the similar stronger damping may 
appear below about 5.5, and 5 T at 0.76 and 0.88 K, 
respectively, where the boundary between the quantum 
vortex slush and the thermal vortex liquid regions is 
located!] This boundary can be seen in the schematic vor- 
tex phase diagram in the low temperature regionta shown 
in the upper left inset of Fig. 5. In the diagram, N is the 
normal, S the vortex solid, VS and TVL the quantum vor- 
tex slush and thermal vortex liquid states, respectively. 
The SdH oscillation amplitude at 0.76 and 0.88 K, how- 
ever, does not show the clear stronger damping down 
to the lowest magnetic field where the SdH oscillations 
can be detected at those temperatures. It means that 
the magnetic fields where the SdH oscillations are ob- 
served are mostly still in the TVL region. Then the mag- 
netic field dependence of the damping seems to follow the 
damping by the thermal fluctuations with temperature. 
This result suggests that the expected stronger damping 
may appear at lower magnetic fields. 

In order to confirm the proposed model on the stronger 
damping in the SdH oscillations in the quantum vortex 
slush state, it is necessary to show the clear relation be- 
tween the stronger damping and the vortex phase di- 
agram. More precise measurements at lower magnetic 
fields in better quality sample are required in future. The 
transport current density dependence of the SdH oscilla- 
tion amplitude is also important to measure in the quan- 
tum vortex slush state because the non-linear behavior 
in the resistivity has been observed. These experiments 
are in progress. 



IV. SUMMARY 

We observed the SdH oscillations on the flux flow re- 
sistivity in the vortex state of the Q2D organic super- 
conductor k-(BEDT-TTF) 2 Cu(NCS) 2 . The additional 
damping of the SdH oscillation amplitude near H C 2 is 
well described by the model on the superconducting (vor- 
tex) fluctuations as well as those observed in the dHvA 
effect. In only the SdH effect, the stronger damping ap- 
pears in the quantum vortex slush region. The stronger 
damping may reflect the perturbation of a phase coher- 
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ence of vortices and quasiparticles in the quantum vortex 
slush state with the short-range order of vortices due to 
moving those by the transport current. 
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